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We study the role of quantum statistics in the self-ordering of ultracold bosons and fermions mov¬ 
ing inside an optical resonator with transverse coherent pumping. For few particles we numerically 
compute the nonequilibrium dynamics of the density matrix towards the self-ordered stationary 
state of the coupled atom-cavity system. We include quantum fluctuations of the particles and the 
cavity held. These huctuations in conjunction with cavity cooling determine the stationary distri¬ 
bution of the particles, which exhibits a transition from a homogeneous to a spatially ordered phase 
with the appearance of a superradiant scattering peak in the cavity output spectrum. At the same 
time the cavity held Q-function changes from a single to a double peaked distribution. While the or¬ 
dering threshold is generally lower for bosons, we conhrm the recently predicted zero pump strength 
threshold for superradiant scattering for fermions when the cavity photon momentum coincides with 
twice the Fermi momentum. 

PACS numbers: 37.30.-hi,37.10.Vz 


I. INTRODUCTION 

Ultracold particles moving in an optical resonator are 
a fast growing research held since the advent of laser 
cooling [1-3]. The establishment of cavity cooling [4] 
extended cooling possibilities to a large class of linearly 
polarisable particles and in the last years cavity optome- 
chanics [5] has experienced a tremendous boost, both the¬ 
oretically and experimentally [3]. Recently, cavity cool¬ 
ing was even extended to the subrecoil regime in the 
domain of ultracold quantum gases where the nonlin¬ 
ear coupled dynamics of held and particles require a full 
quantum description [6-8]. 

Particularly rich physical dynamics appear for 
transversally illuminated particles inside high-Q optical 
resonators (c/. hg. 1), where the cavity-induced backac- 
tion of the scattered held onto the particles induces a 
phase transition to crystalline particle order. First de- 



Figure 1. Polarisable particles moving along the axis of a 
standing-wave optical resonator are coherently pumped in the 
transversal direction. Above a critical threshold the particles 
collectively scatter photons into the cavity while organising 
themselves in the emerging optical potential. 
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scribed semiclassically [9-15] this phenomenon also ap¬ 
pears quantum-mechanically as a phase transition at 
zero temperature [16-22]. Above a critical threshold 
pump laser intensity the atoms form self-sustained regu¬ 
lar Bragg lattices, which maximise the collective scatter¬ 
ing of laser light into the resonator (superradiance). Self¬ 
organisation has been experimentally demonstrated with 
both thermal gases [23, 24] and atomic Bose-Einstein 
condensates (BECs) [25-27]. It bears a close connec¬ 
tion to the generic model of the Dicke superradiant phase 
transition [28-30]. 

Interestingly, the quantum statistical properties of the 
particles have a decisive influence onto the self-ordering 
dynamics. As a particularly striking effect, a resonant 
reduction of the superradiance threshold for fermions al¬ 
most towards zero pump amplitude has been predicted 
when the cavity photon momentum coincides with twice 
the Eermi momentum [31-33]. 

Since we are dealing with an open system with atom 
pumping and photon leakage through the mirrors, fluc¬ 
tuations and nonthermal equilibrium phenomena are es¬ 
sential parts of the full dynamics, as already observed 
experimentally [34, 35]. We consider here the so-called 
dispersive limit at high laser to atom detuning, where the 
particles constitute a dynamic refractive index within the 
resonator and only coherently scatter pump laser light 
into the cavity. The relative phase and magnitude of 
the scattered photons depend on the particle positions, 
while the motion of the latter is governed by the optical 
dipole force exerted by the cavity and pump fields [1]. 
Recent experiments [6, 27, 35] have achieved cavity loss 
rates small enough to be comparable to the particle re¬ 
coil energy scale. Lacking separate timescales, the light 
field then does not adiabatically follow the particle mo¬ 
tion, inducing an even more complex correlated or even 
entangled particle-field dynamics. 

At the semiclassical level and for distinguishable parti- 
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cles the nonthermal properties of the corresponding sta¬ 
tionary equilibrium state and its transient dynamics have 
been studied recently [12-15]. A related study was car¬ 
ried out for bosons at zero temperature with empha¬ 
sis on the effect of cavity losses on the depletion of a 
BEC [36]. In addition, the effect of fermionic statistics 
on the nonequilibrium stationary state has been explored 
with perturbative diagramatic techniques [37]. 

In this work we investigate quantum-statistical effects 
on the nonequilibrium stationary state close to the self¬ 
organisation threshold. Upon including quantum fluc¬ 
tuations and Markovian noise induced by photon losses, 
we compute the density matrix of the joint state of light 
and up to five atoms. Here the particles’ average en¬ 
ergy is not set by a prescribed temperature, instead their 
momentum distribution follows dynamically from the in¬ 
terplay of cavity cooling and fluctuations. In general, a 
nonthermal distribution emerges, which exhibits strong 
atom-field correlations. In principle, even starting from 
a particle ensemble at finite temperature, the combined 
cooling and self-ordering can lead to a quantum degen¬ 
erate ordered final state. 

Since the photons leaking out of the cavity provide a di¬ 
rectly accessible non-destructive probe of the system, we 
characterise the self-organisation transition with the help 
of the Husimi-Q-fnnction, the intensity and the power 
spectrum of the light emitted from the cavity. This pro¬ 
vides an experimentally directly accessible tool, so that 
the whole self-ordering phase transition process can be 
monitored in real time even in a single experimental run 
with minimal backaction [2] . To investigate the full quan¬ 
tum limit of cavity-induced self-organisation and cooling 
we resort to numerical solutions of the master equation 
within a truncated atomic momentum Hilbert space for 
a finite fixed particle number. 


II. EFFECTIVE MODE MODEL 


We consider N ultracold particles within a high-Q op¬ 
tical resonator illuminated by a transverse laser far de¬ 
tuned from any internal resonance (c/. fig. 1). The light 
scattered by the particles into the resonator interferes 
with the pump field and creates an optical lattice po¬ 
tential along the cavity axis. We assume the particles’ 
motion to be restricted to the cavity axis by means of a 
suitable transversal confinement. The coherent time evo¬ 
lution of this joint particle-field system is then governed 
by the second-quantised Hamiltonian [1, 18] 


H = 


J-LI2 


h? d^ 
2m dx^ 


hUoo^a cos‘^{kx) 


-^hr]{a ^ a^) cos{kx) T(x) dx — fiAca’^'a, ( 1 ) 


where k denotes the wave number of the cavity mode 
which also sets the recoil energy scale huj^ = h?k‘^/ ( 2 m), 
m is the mass of a single particle, and L = M2'K/k is 


the unit cell length of the periodic boundary conditions, 
equal to an integer M times the wavelength A = 27r//c. 
The parameter Uq is the light shift per particle, 77 the ef¬ 
fective pump amplitude per atom and a the annihilation 
operator of a cavity photon. The field operators T’t'(x) 
(T(x)) create (destroy) a particle at position x and obey 
bosonic or fermionic commutation or anticommutation 
relations, respectively. The Hamiltonian (1) is expressed 
in a reference frame rotating with the pump laser fre¬ 
quency cjp. Consequently, the detuning Ac := ujp — uJc 
between the pump and the bare cavity resonance fre¬ 
quency explicitly appears in eq. (1). Note that this ef¬ 
fective model describes any kind of small enough linearly 
polarisable particles [ 1 ]. 

As an important feature the Hamiltonian ( 1 ) allows 
for momentum transfers of Ap = Ahk in addition to 
Ap = ±2hk from scattering within the cavity [7, 38]. 
This enables the system to self-order and has also impor¬ 
tant consequences for the cooling behaviour. Note that 
direct particle-particle interactions, which could induce 
arbitrary momentum transfers Ap {e.g. collisions), are 
neglected. 

The second-quantised representation (1) allows us to 
optimally exploit the symmetries of the system and the 
bosonic or fermionic statistics of the quantum particles. 
To this end we expand the particle’s creation and annihi¬ 
lation operators in the complete set of orthonormal mode 
functions 



V Z Y \ sin(nA/ex) 


( 2 ) 


with integer n > 1. The momentum space discretisa¬ 
tion is given by hAk = ^ 21 : jL^ therefore the cavity wave 
vector k is an integer multiple of A/^ {k = MAf^). As 
discussed in [18] the high-energy modes will be excluded 
from the model by introducing a cutoff n < ric. Note that 
the subspaces spanned by odd (sine) and even (cosine) 
parity modes, respectively, are decoupled. For bosons 
initially in a BEC the sine modes remain unpopulated 
and can be excluded from the picture. For a degener¬ 
ate Fermi gas, however, all modes with energies below 
the Fermi energy are filled up, therefore sine and cosine 
modes both have to be taken into account. Nevertheless, 
the dimensionality of the actual problem to be solved can 
be further reduced by exploiting invariant subspaces (c/. 
appendix A). Additionally, we introduce a cutoff for the 
cavity photon number. 

Photons leaking out of the resonator inevitably intro¬ 
duce noise and damping into the system and thus prevent 
a pure Hamiltonian evolution of the system’s Schrodinger 
wave function. Much rather, the system will evolve into 
a stochastic mixture described by the joint particle-field 
density operator governed by the master equation [39] 

p = — [i7, p] + K. (2apa^ — ap — pa^a) , (3) 

with the cavity photon decay rate 2i<i. We neglect spon¬ 
taneous emission of the atoms under the assumption of 
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a large detuning of the laser with respect to any internal 
atomic resonance. 

An expansion of the (bosonic or fermionic) field oper¬ 
ator in terms of the mode functions (2) yields an 

effective mode model, whose explicit Hamiltonian is re¬ 
ported in appendix A. We restrict ourselves to only the 
appreciably occupied modes, to be able to integrate the 
master equation (3) directly. This is possible since the 
occupation of the higher-energy modes remains low due 
to the ongoing cavity cooling effect [9]. 

III. SELF-ORGANISATION IN THE QUANTUM 
REGIME 

Self-organisation can be understood already on a clas¬ 
sical level, i.e. with classical particles and a damped cav¬ 
ity field [9, 11]. At this level, when the threshold pump 
strength for the phase transition is reached, the symme¬ 
try of the system is spontaneously broken as the initially 
homogeneously distributed particles organise and occupy 
one of two possible configurations [1, 9] corresponding to 
even or odd sites of the cavity optical lattice ^ cos^(/cx). 
In the perfectly homogeneous phase, no photons are scat¬ 
tered into the cavity because of destructive interference 
and self-organisation is triggered only by density fluc¬ 
tuations. The occupation of the odd or even potential 
wells, respectively, is associated with one of two oppo¬ 
site phases of the intracavity light field relative to the 
pump laser. The order parameter of the phase transition 
is given by 0 = {cos{kx)) (or equivalently (a)), which is 
zero in the completely homogeneous phase and ±1 in the 
limit that the particles are completely localised at antin¬ 
odes. Quantum-mechanically, the dynamics lead towards 
a self-organised state which is a superposition of the two 
configurations in odd and even wells correlated with the 
corresponding field phase [40]. Initial atom-field entan¬ 
glement decays on the timescale of the photon lifetime 
rendering the stationary state into a mixture of the two 
configurations. 

For both the entangled state and the mixture the mean 
field (a) and the order parameter {cos{kx)) are zero, while 
the photon number (a^a) is finite. Above threshold, the 
stationary state is approximately [16] 

Pss ~ 2 1 ^’ 2 

where ja) is a coherent state with complex field ampli¬ 
tude a and |±) is the state of particles organised in the 
odd or even wells, respectively. With increasing pump 
strength, the components’ coherent states \a) and \—a) 
have less and less overlap. This behaviour can be seen in 
the Q-function [41] of the cavity field, as shown in fig. 2a. 

The classical order parameter {cos{kx))p^^ is always 
zero and thus unsuitable in the quantum case. In order 
to find an alternative measure, we project pss to obtain 
Pa •= (<^|Pss|<^) / tr ((a|pss|<^): with a maximising the Q- 
function obtained from the cavity field associated with 
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Figure 2. a) The order parameter 0 and intracavity photon 
number (a^a), b) the cavity field spectrum |AS'(ct;)| as a function 
of the pump strength, respectively. Insets: The Q-function of 
the photon field splits up into two coherent peaks across the 
threshold. Parameters: A = 2 bosons, Uo/cjr = 0, = 1, 

Ac/cc;r = -1.75. 

Pss- Far above threshold, where (aj — a) ^ 0, we have 
pa ~ 1+) (+1- We then calculate the order parameter 
with this projected particle state, i.e. 

0 := (cos(fca:))p„. (5) 

As long as the Q-function only has one maximum at 
(a = 0, this is exactly zero because of the system’s sym¬ 
metry. The threshold pc? the driving for which 0 
first becomes larger than zero, is therefore defined to be 
the pump strength at which the Q-function splits up and 
has two local maxima. The sudden increase of the mean 
intracavity photon number at the threshold, which is ex¬ 
pected in the thermodynamic limit, is smoothened in the 
few particle regime considered here. This is illustrated in 
fig. 2a. 

In addition to the order parameter, the onset of self¬ 
organisation can also be seen in the cavity field spectrum, 
as observed experimentally with a BEG [26, 42]. To this 
end we calculate S{ijj) = T{{a^{T + t)a{T)))^ where T 
is the Fourier transform with respect to t, and T is a 
time sufficiently large for the system to have reached its 
stationary state [41]. Figure 2b shows |5'(cc;)| for bosons. 
Below threshold, the two peaks correspond to collective 
modes, characterised by light-held oscillations and parti¬ 
cles excited from momentum zero to Hk and vice versa. 
At low pump strength, the peaks are located at ^ Tcjr 
and are almost perfectly sharp (undamped) since the 
collective modes are almost purely atomic excitations. 
With increasing pump strength, apart from a broaden¬ 
ing of the peaks due to mixing of atomic excitations with 
the decaying photonic excitations, we observe a softening 
of the collective modes which become energetically less 
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Figure 3. Time evolution of the momentum population p 
for a) bosons and b) fermions, starting with a quantum de¬ 
generate initial state at temperature zero. Participation of 
the fermion with p = 0 in the photon scattering dynam¬ 
ics is strongly suppressed, as its target states p/(hk) — ±1 
are essentially blocked. Parameters: N = 5, rj/uJn = 2, 
Ac/ct;R = —4, f/o/^R = —0.5, Ak/k = 0.5. 


and less costly. At a critical pump strength these modes 
are shifted to zero energy (u; = 0 in fig. 2b). This pro¬ 
cess, accompanied by collective scattering of photons into 
the cavity at the laser frequency, is called superradiant 
because of the quadratic scaling of the photon number 
with N. The presence of such a dominant coherent peak 
at cj = 0 is thus an additional characterisation of self¬ 
organisation, as already observed semi-classically [15]. 
Within our model this superradiant coherent peak has 
zero width since we do not include a finite linewidth of 
the pump laser. 

For fermions, the behaviour of the order parameter 
is qualitatively similar to the bosonic case illustrated in 
fig. 2a. The spectrum, however, shows qualitative differ¬ 
ences which will be discussed below. 


IV. QUANTUM STATISTICS AND 
SELF-ORGANISATION 

Self-organisation and superradiance have been demon¬ 
strated with both a thermal gas [23, 24] in the classical 
regime as well as for a BEC close to T = 0 [25, 27]. 
Self-ordering of course requires excitations out of the 
p = 0 condensed state to higher momenta with incre¬ 
ments (c/. fig. 3a). The qualitative behaviour of 

BEC self-organisation resembles its thermal counterpart, 
with the recoil energy scale huj^i playing the role of the 
temperature scale /cpT [43]. 

In contrast to bosons, fermionic statistics allows for 
a qualitatively different self-ordering behaviour [31-33]. 
In fact, a Eermi gas introduces a new internal energy 
scale, the Eermi energy ep, and indeed its ratio to the re¬ 
coil energy hujji dramatically affects the behaviour of the 
fermionic self-organisation. As illustrated in fig. 3b, due 
to the Pauli pressure fermions occupy higher momentum 
states even at T = 0 including cosine and sine modes. 



I ' ■ ■ I ■ ' ■ I H N I ' ' ■ I ■ ■ ' l -l LJ 0.0 

-2 -1 0 - 2-1 0 
(Ac - AT/o/2)/u;r (Ac - NUo/2)/uj^ 


Figure 4. The phase diagram of two and three bosons and 
fermions, respectively. The threshold is compared with the 
analytic results of [12] in the semiclassical limit (dashed) 
and [43] for bosons (dash-dotted), as well as [37] for fermions 
(solid, kF/k = 1). For bosons, the sudden decrease of the 
order parameter with increasing pump strength is an artefact 
of the photon number cutoff. Parameters: Uq/luk = —0.5, 
I^/UJR = 1 . 


This prevents a direct mapping of the system’s Hamilto¬ 
nian onto an effective Dicke model (demonstrated for a 
BEC [44]) and the understanding of self-organisation in 
terms of a simple superradiant transition is not possible, 
except for the limiting case of an infinitely large Eermi 
momentum relative to the cavity photon momentum [45] . 

The central effects of Eermi statistics on self-ordering 
so far have been evaluated based on the assumption of 
a thermal equilibrium for the particles and neglecting 
the role of quantum fluctuations of the cavity field [31- 
33]. Here we test these predictions in comprehensive 
nonequilibrium simulations for a small number of par¬ 
ticles including light-matter entanglement and quantum 
fluctuations of the cavity field. Due to Pauli blocking we 
can expect that for fermions not all particles will equally 
participate in the self-organisation phenomena. Whilst 
for classical particles and bosons scattering of a photon 
is always allowed and thus the critical pump strength 
scales as pc oc 1/\/N [10, 12, 21, 43], only fermions cou¬ 
pled to an unoccupied final state can scatter (c/. fig. 3b). 
This fraction strongly depends on the Eermi energy and 
hence scales differently with N. Eor a sharp Eermi sur¬ 
face, Pc can be estimated analytically. In one spatial di¬ 
mension one obtains ep oc {N/L)‘^, which in turn means 
Pc oc ^y\n(N)JN for /c ^ /cp or even pc oc a/TV for 
k k^ [31]. 

Eigure 4 depicts the value of the order parameter (5) 
as a function of the (dispersively shifted) cavity detuning 
and the pump amplitude for two and three bosons and 
fermions, respectively. In this low-particle regime the 
notion of a sharp Eermi surface becomes questionable. 
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Figure 5. a) Cavity field spectrum and b) initial conditions with lowest-energy particle excitations for bosons (I) and fermions 
(II-IV). For fermions the value Ak/k = 1,0.5,0.25 is varied. For bosons only the case Ak/k = 1 is shown, as the dynamics 
are unchanged in the other cases. Note that in case II the peaks in the spectrum at dicjR correspond to transitions which 
are forbidden in the initial state but become possible during relaxation towards the stationary state. The frequency of the 
transition indicated by the arrow in II b) lies outside the depicted spectral range. The resonance condition k = 2/cf is fulfilled 
in IV. Here the central peak in the spectrum at the laser frequency (lu = 0) is dominant even for the smallest pump strengths. 
Parameters: V = 4, Uq/cj-r = —1/16, k/ujr = 1/8, Ac/ujr = —1/2. 


Nevertheless, we see qualitative agreement of our numer¬ 
ical data with the scaling discussed above: for = 3 
fermions the rescaled yfNrjc is larger than for A/" = 2, 
while for bosons it does not change. 

While the Fermi blockade causes an unfavourable 
threshold scaling with the number of particles, the ex¬ 
istence of a Fermi surface can also give rise to a resonant 
reduction of the threshold if k = 2/cf. This “umklapp 
superradiance” [31] involves particles scattered from one 
end of the Fermi surface to the other (ie. from —k^ to 
k^) by the two-photon transition with momentum trans¬ 
fer hk = 2hk-F. This process has essentially no energy 
cost and thus can drastically lower the self-organisation 
threshold. Here the ordered state is very close to the 
Fermi surface and has essentially the same kinetic energy 
as the homogeneous ground state. Note that this avail¬ 
ability of some already excited particles can also explain 
the threshold reduction found in [21] for bosons at small 
temperature compared to the zero temperature limit. 

In our simulations with a fixed number of fermions pre¬ 
pared initially at T = 0, the resonance condition k = 2kY 
can be reached by different means: either by adjusting 
the cavity wavenumber k directly or by altering the par¬ 
ticle density L/N and thus the Fermi momentum. Here 
we choose to keep k fixed and vary the trap volume L 
which the fermions can occupy by altering Ak = 27r/I/. 
This has the advantage that the recoil frequency cjr stays 


the same for all parameters. 

In order to investigate the role of umklapp superradi¬ 
ance, we prepare the initial state in the resonant condi¬ 
tion, perform the full quantum dissipative evolution and 
compute the order parameter and intracavity spectrum 
once the stationary state is reached. As mentioned in 
the beginning, our simulations include the relaxation of 
the particles into a (generally) nonequilibrium station¬ 
ary state different from the initial one. In particular, 
redistribution of energy and momentum will alter the 
Fermi surface, thereby modifying the resonance condi¬ 
tion. Whether and how the umklapp superradiance man¬ 
ifests itself in this stationary state is a non-trivial ques¬ 
tion. The answer is illustrated in fig. 5. Here we compare 
bosons in a BEC with zero-temperature fermions at three 
different densities corresponding to the three Fermi mo¬ 
menta k^ = 2k {Ak = k), k^ = k {Ak = 0.5k) and 
2k-F = k {Ak = 0.25k). In the last case the resonance 
condition is fulfilled. Initially, with no cavity photons 
present, the particles move freely with a quadratic dis¬ 
persion relation. As already mentioned, the cavity only 
supports momentum transfer to and from the particles 
in units hk. This corresponds to steps of 1 on the x-axes 
of fig. 5b, where arrows indicate the particle excitations 
with lowest possible energy cost. These excitations are 
visible in the spectra shown in 5a as peaks at the corre¬ 
sponding frequencies. 






















6 



VNtj/uj-ei 


Figure 6. a) Order parameter and b) photon number for 
fermions (solid) and bosons (dashed) for the resonance case 
k = 2/cf- The threshold for fermions is lower compared to 
bosons. Insets: The Q-function of the photon field at a pump 
strength which lies above threshold for fermions and below 
threshold for bosons, respectively. Same parameters as in 
fig. 5. 


tallisation of quantum particles in cavity fields even in 
this nonequilibrium context. This can be nicely moni¬ 
tored in real time with minimal perturbations by observ¬ 
ing the spectrum of the scattered light. 
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Appendix A: Mode model 


The effect of resonant umklapp scattering is clearly 
visible in the spectra by comparing the initially reso¬ 
nant case (fig. 5a IV) with the other two fermionic non¬ 
resonant cases as well as with the bosonic case. The 
appearance of a coherent peak at zero frequency in the 
resonant case signals the onset of superradiance even at 
almost vanishing pump strengths. Eor all the other cases 
this onset of a coherent scattering peak is located at 
higher pump strengths. Such a low superradiant thresh¬ 
old is possible because of (almost) resonant particle exci¬ 
tations with momentum hk, creating the correct density 
modulation period for coherent Bragg scattering into the 
cavity. 

Even though the umklapp resonant condition can lower 
the threshold for the appearance of the superradiant 
peak much below the bosonic threshold, the number of 
photons in the cavity and the amplitude of the order 
parameter (5) can still be lower than in the bosonic 
case, due to the unfavourable scaling with N illustrated 
above. Eor instance, the threshold pump strength scales 
as rjc ^ \/N/ ln(A) at k 2/cf if we assume a sharp 
Eermi surface. Therefore the effect of umklapp super¬ 
radiance is much less visible in the photon number or 
the order parameter as compared to the spectrum. Still, 
for certain parameter regimes and in the resonant case a 
lower self-organisation threshold for fermions can be seen 
even in these quantities, as shown in fig. 6. 


Starting from the Hamiltonian (1), we use the mode 
model from [18], extended to include the sine modes and 
generalized for an arbitrary momentum space discretisa¬ 
tion A/c. 

We expand the particles’ field operator ^(x) in the 
complete set of orthonormal mode functions 







cos(nA/cx), 





sin(nA/cx), 


(Al) 


with integer n > 1. The modes are labeled by the multi¬ 
index j = (j, Jtt), which contains the mode number j and 
the parity = 1 for the cosine and = — 1 for the sine 
modes. This yields the decomposition 


^i.x) (A2) 

3d<nc 


where the Cj follow either bosonic or fermionic (anti-) 
commutation rules and the Cj^-i operators are formally 
set to zero for bosons. 

With c = (co,i,.. ■ as column 

vector and = (4,i’• • • ’ 4e,i’4-i> • • • > -i) as row 

vector, we can write the Hamiltonian (1) as 


V. CONCLUSIONS 

Self-ordering of ultracold particles in an optical res¬ 
onator field survives as a stationary nonthermal equilib¬ 
rium state in the long-time limit even if cavity losses with 
their inherent field quantum fluctuations are taken into 
account. The quantum statistics of the particles turns 
out to have a decisive influence on the light-induced crys- 


H = —hAca^a^huj^i Kc^ ^hr]{a^a^) 

+ (ct (J + C2fe/A J a) , (A3) 

where cjr = Kk^ j(2m) is the recoil frequency and I is 
the identity matrix. The matrices K and Cm with inte¬ 
ger M, producing the kinetic energy term and the terms 
proportional to cos{kx) and cos(2/cx), respectively, are 
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defined as 



2" 




(A4) 


where a = 1/v^ and b = 1/2. For Cm^ the two blocks 
belong to the cosine and sine modes, respectively. The 
two populated secondary diagonals are separated by M 
modes. Additionally, modes j, k for which j = M are 
coupled by the entries b for cosine modes and —b for sine 
modes, respectively. For reference, the quadratic forms 
appearing in the mode-expanded Hamiltonian (A3) read 


Cm^ — ^ 0,1 

[Sj,k+M + Sj+M,k + jTrSj+k,M] C^jCk- (A5) 


^ j,k 


For our simulations we use the usual basis of Fock states 

3,15 • • • : ^nc,l 5 ^1,-15 • • • 5 ^nc, —1) f 5 ('^^) 

^ nn=N 


where rij = 0,..., for bosons and rij = 0,1 for 
fermions. 


As already noted and as can be seen from the cou¬ 
pling matrices (A4), sine and cosine modes are decou¬ 
pled. However, for M := k/Aj^ > 1 one can identify 
further subgroups of modes which are decoupled from all 
other groups. We will now specify these groups and use 
them to identify invariant subspaces of the Hilbert space 
H spanned by the basis states (A6). 

Let us group all available modes into s sets of modes, 
where s := [M/2J + 1 for bosons and s := 2[M/2j + 1 
for fermions: 

9i ■= {j\jn = ln,j mod M = l\J j mod M = M -l] , 

(A7) 

with = 1 = ±1) for bosons (fermions) and 

I < [M/2J. The significance of these sets is that two 
modes are coupled through the dynamics induced by the 
Hamiltonian (A3) if and only if they belong to the same 
set. As a consequence, the number operator correspond¬ 
ing to each group, i.e. 


ni-.= Y^ c\cj (A8) 

3^91 

commutes with H and {hi) is constant over time. This 
enables us to define invariant subspaces of H: two basis 
states of (A6) belong to the same invariant subspace if 
and only if they produce the same s-tuple {{hi))i. We 
can exclude any basis state from our model if it belongs 
to a subspace which is not populated by the initial state. 

For example, for a degenerate Fermi gas with A" = 4, 
k/Ak = 4 and ric = 6, the basis (A6) consists of 715 
states. By virtue of the above considerations we only 
have to take 72 basis states into account. 
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